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Abstract 

With sintered needles aligned and a magnetic field applied transverse to its easy axis, the rare- 
earth ferromagnet Nd 2 Fei 4 B becomes a room-temperature realization of the Random Field Ising 
Model. The transverse field tunes the pinning potential of the magnetic domains in a contin¬ 
uous fashion. We study the magnetic domain reversal and avalanche dynamics between liquid 
helium and room temperatures at a series of transverse fields using a Barkhausen noise technique. 
The avalanche size and energy distributions follow power-law behavior with a cutoff dependent on 
the pinning strength dialed in by the transverse field, consistent with theoretical predictions for 
Barkhausen avalanches in disordered materials. A scaling analysis reveals two regimes of behavior: 
one at low temperature and high transverse field, where the dynamics are governed by the ran¬ 
domness, and the second at high temperature and low transverse field where thermal fluctuations 
dominate the dynamics. 

BAGS numbers: 75.60.Ej,75.50.Lk,75.60.Jk,75.10.Hk 
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I. INTRODUCTION 


The response of a ferromagnet to a time varying magnetic field has been stndied exten¬ 
sively via experiment, nnmerical modeling, and theory.^^® The magnetization, rather than 
evolving continuonsly as the field changes, often increments in discrete steps as individnal 
domains or sets of domains rapidly switch direction. These steps were hrst observed as volt¬ 
age pnlses in indnctive picknp coils wrapped aronnd the magnet, picking np the change in 
magnetization with time, dM/dt. The crackling noise arising from the snperposition of many 
snch pnlses is known as Barkhansen noise.^ Detailed measnrements of the spectrnm of this 
noise provide a sensitive probe of the internal energetics of domain reversal,^ with nses both 
in industrial non-destructive evaluation of magnetic materiald^ and as a valuable technique 
for understanding domain dynamics for model magnetic materials and Hamiltonians.^ 

One important model system amenable in principle to such measurements is the Ran¬ 
dom Field Ising Model (RFIM), which is considered to be a general model for disordered 
systems.^ In the RFIM, the standard Ising Hamiltonian is modihed by adding an applied 
field along the Ising axis whose value varies randomly from site to site: 

« = + ( 1 ) 

ij i 

where is a site-random field and J is the interspin interaction. The original experimental 
studies of the RFIM followed a proposal by Imry and Ma,^ in which a site-diluted Ising 
antiferromagnet in a large static magnetic held forms a realization of the RFIM Hamil¬ 
tonian of Eq. 1 . While this approach proved fruitful for studying quantities such as the 
thermodynamic critical exponent^^^^ and correlation lengths,^ the lack of a net long wave¬ 
length moment in antiferromagnets limits the potential probes and hence the set of physical 
questions that can be addressed. 

Within the last decade, ferromagnetic realizations of the RFIM have been discovered. 
They arise when diluted, dipole-coupled Ising materials are placed in a static magnetic 
held applied transverse to the Ising axis. The combination of the oh-diagonal elements 
of the dipole-dipole potential and the symmetry breaking of the disorder lets the uniform 
transverse held generate an ehective site-random longitudinal held. The initial realizations, 
hrst in the rare-earth huoride LiHoa;Yi_3;FiiSHIll and subsequently in the molecular magnet 
Mni2-acj ^^ * ^° * used dipole-coupled single spins. Both materials exhibit magnetic ordering at 
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cryogenic temperatures due to the low energies of the dipolar coupling. Recently, Tomarken 
et al. demonstrated^ that a sintered block of the rare-earth ferromagnet Nd 2 Fei 4 B acts as a 
realization of the RFIM at room temperature. In Nd 2 Fei 4 B, individual spins are exchange- 
coupled and form elongated domains due to the inherent crystalline anisotropyj ^^F3 | These 
domains then interact via a dipolar coupling. The extended nature of the dipoles gives 
them much larger moments than individual spins would have, with a correspondingly higher 
energy scale (e.g. the Curie temperature for sintered Nd 2 Fei 4 B is 585 hPl). The random 
packing of the grains from the sintering process yields the necessary translational disorder 
for RFIM behavior. 

We exploit the ready availability and tunability of the RFIM state in Nd 2 Fei 4 B to record 
the hrst experimental measurements of Barkhausen noise in a random-held magnet. As 
the transverse held deepens the pinning wells, the power law distributions in avalanche size 
and energy are truncated, consistent with simulation results. Critical exponents are close 
to mean held values in the low temperature, high transverse held regime, where we also 
discover extended oscillatory events in the magnetization caused by initial overshoot of the 
hnal equilibrium state by large-domain reversals. 


II. METHODS 

We measured the static magnetization of 2 mm diameter by 10 mm length cylinders of 
commercially-sintered Nd 2 Fei 4 B via Hall magnetometry, using a pair of passive GaAs Hall 
sensors (Toshiba THS118) in a gradiometric conhguration to remove the contribution of the 
applied held to the Hall signal. The easy axes of the samples were oriented parallel to the 
cylindrical axis. We measured the domain dynamics via Barkhausen noise measurements. 
A 320 turn pickup coil was wrapped around the same Nd 2 Fei 4 B cylinders and an identical- 
geometry, empty coil wired in opposition to cancel the dH/dt from the applied held. The 
signal was amplihed by 10,000 with a Stanford Research SR560 low-noise voltage pream- 
pliher, followed by a 120 kHz Krohn-Hite 3988 low-pass hlter to eliminate aliasing ehects, 
and hnally a National Instruments USB-6211 250 kHz/16 bit digitizer. Longitudinal and 
transverse magnetic helds were supplied by a dual axis 5 T/2T superconducting magnet. 
Hysteresis loops were typically obtained with the transverse held hxed and the longitudinal 
held ramped at a constant 0.4 T/min between ± 4.5 T. Given the large held scales involved. 
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it was not necessary to screen out the Earth’s magnetic held. 

For each temperature/transverse held point investigated, a series of hysteresis loops was 
measured to obtain a minimum of 10® events; depending on temperature and held, this 
required 25 to 38 loops. For measurements at T = 4.2 K, helium exchange gas was used to 
place the sample in thermal contact with the liquid helium bath; at T = 150 and 300 K, the 
sample was in vacuum and a PID loop was used to stabilize the temperature. After acqui¬ 
sition, a software 1.1 kHz high-pass hlter was applied to the data to remove static ohsets, 
external power-line noise, and other spurious low-frequency ehects. Individual events were 
identihed by looking for voltages that deviated from the mean by more than 3.5 standard 
deviations of the time-independent background-volt age noise of the measurement chain. 


III. RESULTS 

We plot in Fig. [^a) the overall magnetic behavior of Nd 2 Fei 4 B at T = 150 K, well 
below the 585 K Curie temperature.^^ Applying a constant 6 kOe transverse held results in 
several signihcant changes to the hysteresis loop. First, there is a decrease of approximately 
5% in the saturation magnetization, due to the hnite value of the intrinsic anisotropy and 
consequent tilting of spins away from the Ising (easy) axis. Second, the shape of the loop 
changes when the transverse held is applied. As discussed in detail in Ref. EU a transverse 
held decreases the enclosed area of the hysteresis loop for longitudinal helds below 15 kOe, 
while for larger longitudinal helds, applying a transverse held increases the loop area. This 
broadening is a signature of enhanced pinning in the high-held regime where the typical 
domain size is large due to coarsening, also substantiated by studies of Return-Point Memory 
trajectories.^^^^ 

In order to study the domain reversal dynamics in the weak and strong pinning regimes, 
we performed Barkhausen measurements across full hysteresis loops for requisite transverse 
helds and temperatures. A schematic of the energy landscape is shown in Fig. [^b), where 
the combination of intrinsic and extrinsic pinning with the constantly ramping longitudinal 
held produces a tilted washboard potential for the magnetization state of the system. In the 
weak-pinning limit (top row), the relatively shallow minima allow the system to progress 
monotonically through the set of minima, resulting in a series of single-signed pulses in 
dM/df and hence single-signed spikes in the induced voltage. This is the behavior typically 
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FIG. 1. Overview of random-field effects in Nd 2 Fei 4 B. (a) Hysteresis loops at T = 150 K, with 
magnetization measured parallel to the Ising (easy) axis in static 0 and 6 kOe transverse magnetic 
fields. Application of a transverse field changes the shape of the loop, narrowing it at low longi- 
tndinal field and broadening it at high longitudinal field. The latter is an indication of enhanced 
pinning. Boxes at low and high longitudinal field show the ranges for the Barkhausen data plotted 
in Figs. and respectively, (b) Schematic showing the evolution of the free energy landscape 
when a transverse field is applied. The transverse-field-induced random field tunes the system from 
the weak to the strong disorder limit, resulting in deeper local minima in the energy landscape, 
where the magnetization can oscillate about these minima. 

observed in Barkhausen experiments on magnetically soft and amorphous materials such as 
ferroglasses.^^HEl gy contrast, increasing the depth of a number of the pinning wells (bottom 
row) can result in qualitatively different behavior. It becomes possible for a large reversing 
domain to overshoot its ultimate equilibrium point and oscillate about the center of the 
potential well. This results in transitory oscillatory behavior in the magnetization as the 
system comes to equilibrium, and hence oscillatory behavior in the pickup coil voltage. 
Behavior of this kind was recently reportecP in a steel alloy in which uniaxial stress was 
used to enhance the strength of the pinning. 

The salient features of our experiments are captured in Figs. and These plots 
present raw time series of the voltage induced in a pickup coil surrounding a Nd 2 Fei 4 B 
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FIG. 2. Time series of induced voltage on a pickup coil wrapped around a Nd 2 Fei 4 B cylinder 
at T = 150 K, low longitudinal field (increasing continuously at 0.4 T/min), and a series of 
transverse fields, (a) At low longitudinal fields, the absence of domain coarsening gives time series 
dominated by single-signed pulses, corresponding to monotonic changes in the magnetization, (b), 
(c) Successive enlargements of one reversal event at 6 kOe transverse field, showing the pulse 
structure over a duration of approximately 8 ^s. 

cylinder for two segments of the hysteresis loop at T = 150 K (Fig. [^a)) in three hxed 
transverse helds. At the low longitndinal helds of Fig. the system consists of an ensemble 
of small domains that reverse monotonically as the longitudinal held is ramped, following the 
behavior sketched in the top row of Fig. [^b). The rate of occurrence of visible Barkhausen 
spikes drops as the strength of the disorder increases with increasing transverse held. This 
suggests that enhancing the disorder-induced pinning in this regime suppresses avalanches 
consisting of multiple domains reversing in quick succession to produce an observable voltage 
spike and forces the system towards a magnetization change dominated by small domains 
reversing individually, causing voltage spikes which are at or below the noise hoor of our 
measurement chain.®! As shown in the successive enlargements of Figs. ib) and (c), the 
observable events typically last less than 10 microseconds in duration. The fast dynamics 
are a result of the strong pinning in general in Nd 2 Fei 4 B; in the ferroglasses and other soft 
materials more commonly investigated using Barkhausen techniques, the reversal timescales 
are typically 2 to 3 orders of magnitude longer.!2®2Sl 

At the high longitudinal helds of Fig. we observe different signatures in the Barkhausen 
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FIG. 3. Time series of induced voltage on a pickup coil wrapped around a Nd 2 Fei 4 B cylinder 
at T = 150 K, large longitudinal field (increasing continuously at 0.4 T/min), and a series of 
transverse fields. In contrast to the behavior at low longitudinal fields (cf. Fig. [^, the deep 
potential energy wells created by the combination of domain coarsening and the random-field- 
induced pinning enable magnetization oscillations, (a) Time series for a range of transverse fields. 
In the absence of a transverse field, the system is in the weak disorder limit and evolves in a nearly 
continuous fashion. As the transverse field is increased, the strength of the disorder increases, 
giving a steadily increasing number of discrete events, (b), (c) Successive enlargements of a single 
event at a longitudinal field of 3.2 T. The overall duration of the oscillatory event is approximately 
30 ms, with an oscillation period of 44-48 fis (11-12 sampling points). 


noise spectra. Domain coarsening effects result in large domains so that even the reversal 
of an individual domain can produce a detectable signal. As the strength of the pinning is 
again increased by ramping the transverse held, we now see an increasing frequency of events 
where the magnetization oscillates repeatedly, corresponding to the physical picture at the 
bottom of Fig. [^b). The typical oscillation period is of order a few tens of microseconds, 
with the ringdown period of the entire event lasting for a few tens of milliseconds (Figs, [^b) 
and (c)), indicating a relatively low degree of anharmonicity in the local pinning potential. 
We note that even in the low temperature limit the application of a transverse held is still 
required to induce magnetization oscillations. 
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IV. DISCUSSION 


The dynamics of the monotonic switching events of the type plotted in Fig. can be 
investigated by examining probability distribntion histograms of different moments of the 
events and stndying the evolution of the distributions as a function of temperature and 
transverse held. These probability distributions are shown in Fig. |^for the integrated area 
J(v- r„) dt (Fig. 08)) and the total event energy E = f(V — Vq)^ dt (Fig. |^b)), 
where V is the amplihed time-dependent voltage and Vq is the measured noise hoor of the 
ampliher chain. S is proportional to the total change in magnetization associated with 
a given reversal event, whereas E measures the energy dissipated over the course of the 
event 

We observe power-law behavior for both S and E. As a function of temperature and 
transverse held, the critical exponents fall into two groups. At T = 4.2 K and = 3 and 
6 kOe, and at T = 150 K and Ht = 6 kOe, the exponents are 2.3 ± 0.2 and 1.9 ± 0.2 
for the distributions of S and E, respectively. The exponents for S and E are comparable 
to the values 9/4 and 11/6 predicted from mean-held theorj^^^^, as would be expected for 
a long-range interaction such as dipole-dipole coupling. Moreover, the energy distribution 
demonstrates a power-law relationship to the area (Fig. |^b) inset) with exponent 1.64±0.1, 
not far from the value of 1.5 predicted from mean-held theory.!^ For all transverse helds at 
T = 300 K, as well as the low held behavior at lower temperatures, we hnd a diherent set 
of critical exponents: 3.5 ± 0.4 for S and 2.4 ± 0.3 for E. This bimodal behavior suggests 
that we can divide the system into two regimes, the hrst where the disorder induced by the 
transverse held dominates over thermal ehects and the second where the thermal behavior 
dominates. In the disorder-dominated regime, we see behavior broadly consistent with the 
mean-held expectation. 

We compare our measurements with the results of numerical simulations of the RFIM by 
examining a scaling collapse of the probability distributions of the event sizes for diherent 
temperatures and transverse helds. Following the approach of Ref. |3Sl we attempt to scale 
the probability distributions using the functional form D^rS"") = R), 

where a is the exponent describing the cutoh threshold, r + af36 parameterizes the slope of 
the probability distribution integrated around a complete hysteresis loop above the cutoh, 
and r = {R — R^) /R is a reduced disorder.®! We plot in Fig. |^a) contours for the scaling 
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FIG. 4. Probability distributions of avalanche energies and sizes for T = 4.2, 150, and 300 K, and 
transverse fields Ht = 0, 3, and 6 kOe, in the low longitudinal field regime plotted in Figurej^ Solid 
lines are fits showing the extent of pure power-law behavior; dashed lines are extrapolations of those 
fits, showing deviations from the power laws. Traces are offset vertically for clarity, (a) Probability 
distribution of integrated event area S. The exponents top to bottom are 3.1(0.4), 2.2(0.2), 2.5(0.3), 
3.5(0.4), 3.6(0.5), 2.5(0.3), 3.5(0.4), 3.7(0.8), and 3.3(0.6). (b) Probability distribution of event 
energies E. The exponents top to bottom are 2.2(0.3), 2.0(0.2), 1.9(0.2), 2.3(0.2), 2.8(0.4), 1.9(0.2), 
2.4(0.2), 2.3(0.7), and 2.5(0.5). Inset: avalanche size versus energy at 77* = 6 kOe and T = 4.2 K. 
Solid line is a power-law fit with exponent 1.64 ± 0.1. 

of all nine temperature/transverse field values as a function of the two scaling exponents. 
We see that the best-fit value for the cutoff threshold exponent a is 2/3 for all fits, with only 
a weak dependence on temperature and transverse field. This value exceeds the mean-field 
prediction of 1/2, but lies closer to it than the value of 1/4 calculated for nearest-neighbor 
interactions in d=3.^ By contrast, we find that the best-fit values for r -|- a/3S fall into two 
classes depending on temperature and transverse field. 

In the disorder-dominated regime (transverse field large compared to temperature), the 
probability distributions can be scaled on top of each other with little scatter (Figs, [^b) 
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FIG. 5. Scaling behavior in the strong disorder limit, (a) contours for the avalanche area scaling 
parameters a and r + cr/JJ. White circles mark the locations of the best-fit minima, with deviations 
ranging from small (blue) to large (red). Two distinct classes of behavior are observed. At large 
Ht and small T, where domain pinning is strong, there is a narrow minimum at cr = 0.67 ± 0.02 
and T -|- a 136 = 2.6 ± 0.25. At smaller helds and/or larger temperatures, where pinning cannot 
compete with thermal fluctuations, qualitatively different scaling behavior is observed. There is a 
well-defined minimum at r -|- (t/ 3<5 = 3.75 ± 0.4 and a broader dependence on a with the absolute 
minimum remaining at 0.67 ± 0.04. (b), (c) Scaling collapses for the size and energy distributions 
in the strong disorder limit (T = 4.2 K, Ht = 3 and 6 kOe, and T = 150 K, Ht = 6 kOe). The solid 
lines are a phenomenological function from Ref. |33l The scaling parameters for the size distribution 
(b) were a = 0.67 and r -|- a/36 = 2.53, 2.47, 2.70 from top to bottom and for the energy (c) were 
aE = 0.33 and {T + a^6)E = 2.2, 2.1, 2.1. 


and (c)). In this regime, at T = 4.2 K and Ht = 3 and 6 kOe, as well as T = 150 K and 
Ht = 6 kOe, the area probability distribntions have best-£t valnes of r -|- a/36 = 2.6 ± 0.25, 
where as reported above a = 0.67 ± 0.02. This valne of r -|- a/36 compares to the mean- 
held expectation of 2.25 and the nearest-neighbor valne of 2.03.^ A similar scaling collapse 
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analysis on the energy distribntion for the same set of temperatnres and transverse helds 
yields aE = 0.33 ± 0.04 and (r + (7(36)e = 2.1 ± 0.2, compared to mean-held predictions of 
1/3 and 11/6 respectively.*^ The probability distribntions for this disorder-dominated regime 
can be scaled on top of each other with little scatter (Figs, [^b) and (c)). By contrast, for all 
transverse helds at T = 300 K, along with 0 and 3 kOe at T = 150 K and 0 kOe at T = 4.2 K, 
r -f a/36 = 3.7 ± 0.4. In this thermally-dominated limit, the observed best-ht exponents are 
well away from the predictions for both mean-held and nearest-neighbor interactions. The 
scaling is of marginal qnality, snggesting that this regime is best nnderstood as a qnalitatively 
diherent regime from the randomness-dominated regime and that a strong-disorder model 
is not appropriate. 

Importantly, the power-law behavior does not extend to the limits of the probability 
distribntion. As can be seen from the dashed lines in Fig. which are extrapolations of the 
power-law hts, the measnred histograms fall oh from the power law at large S and E. This 
cntoh in the size and energy distribntions of avalanches is a signatnre of a disorder-dominated 
system.*2^*23] Pqj. RFIM in particnlar, the cntoh is expected to follow 5max ~ (7? — 
where S'max is the size of the largest avalanche for which power-law behavior is expected, 
R is a measnre of the strength of the disorder, Rc is a critical level of disorder, and the 
exponent cr is 1/2 in the mean-held limit 

We plot in Fig. |^the measnred deviations from power-law behavior in the energy spec- 
trnm of the avalanches at T = 300 K for a series of transverse helds. Given that we are 
testing the ehects of disorder, we inclnde in this analysis avalanches from both low (Fig. 
and high (Fig. longitndinal helds; qnalitatively similar conclnsions can be drawn from 
considering only the single-signed Barkhansen signatnres of Figs. and In the absence of 
a transverse held, the disorder ehects in Nd 2 Fei 4 B are relatively weak, and the power-law 
behavior extends ont to the highest avalanche energies. When a transverse held is applied, 
the strengthening of the pinning associated with the ehective random held yields a larger 
ehective disorder, resnlting in a cntoh from the power-law at systematically lower avalanche 
energies. This experimental resnlt dovetails nicely with the predictions from the nnmerical 
simnlations of the ehects of increasing disorder.^EHS] 

Finally, we consider separately the magnetization reversal events prevalent at large lon¬ 
gitndinal helds. Instead of a monotonic reversal of a domain or series of domains resnlting 
in a single-signed peak in the indnced voltage, in this regime we see events consisting of 
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FIG. 6. Departures from power-law behavior in the avalanche energy distribution, (a) T = 300 K 
avalanche energy distribution for a series of transverse helds fitted to a power law; traces are offset 
vertically for clarity, (b) Fractional deviation from power-law behavior, showing measured proba¬ 
bility distribution divided by the individual fitted power laws. Transverse-field-induced randomness 
causes systematically more profound deviations from the power-law form. 


a well-defined sinusoidal oscillation which rings down over the course of 1000 or more pe¬ 
riods (Fig. 1^. As the induced voltage is proportional to dM/dt, this is indicative of an 
oscillation in the magnetization itself. Such an oscillation requires a well-dehned harmonic 
minimum in the free energy landscape, along with a signihcant degree of overshoot on the 
initial transition into that minimum. The hrst requirement suggests that the prevalence 
of such oscillatory events will be enhanced with increasing transverse held and hence the 
increasing strength of random-held-induced pinning. The large amount of inertia associated 
with the second requirement limits this effect to large domains reversing as a single unit, i.e. 
primarily in the large longitudinal held regime where coarsening has increased the typical 
domain size. 

In order to assess whether there is a single energy scale associated with these oscillatory 
events, we show in Fig. probability distributions for the oscillation period for T = 4.2 
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FIG. 7. Distribution of periods for oscillatory magnetization events at T = 4.2 K and a range of 
transverse fields. Increasing the transverse field increases the rate of occurrence of the oscillation 
events, while the characteristic oscillation period remains constant. Inset; Normalized distribution 
of oscillation periods for T = 4.2 and 150 K and Ht = 3 and 6 kOe, demonstrating that the shape 
of the distribution is essentially independent of both temperature and transverse held. Smooth 
curve is a Lorentzian ht simultaneously to all four temperature/held sets. 


K and a series of transverse fields. There are two characteristic periods for the oscillation 
events, with a small nnmber of events exhibiting a 16 /is (4 pixel) period and the majority of 
the events showing periods ranging from 30 to above 100 /is, with the predominant period 
at 44 /is. With the application of a transverse field and stronger pinning, there is marked 
increase in the occnrrence rate. This occurrence rate is the primary avenue for temperature 
and transverse-held dependence (Fig. inset). The peak period and shape of the distri¬ 
bution curve are largely independent of the external variables, suggesting that the ratio of 
typical domain size in the coarsening limit to the depth of the pinning well is essentially 
constant and both quantities scale similarly as the disorder strength is tuned upwards via 
the transverse held. A small number of short-period oscillation events are also observed at 
low longitudinal helds. Unlike the longer-period events in the high-longitudinal held regime, 
these overshoots have a typical lifespan of approximately 2 periods, presumably due to small 
domains oscillating in much softer harmonic wells. When a transverse held is applied and 
the strength of the pinning wells is increased, these short oscillation modes are progressively 
suppressed and monotonic reversal events dominate the probability distribution. 
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V. CONCLUSIONS 


While the avalanche dynamics of the Random Field Ising Model have been stndied ex¬ 
tensively via theoretical modeling and numerical simulations, the difficulty of realizing the 
RFIM in bulk ferromagnets has limited the experimental possibilities. The aim of this work 
has been to begin to £11 that gap. We have studied the reversal dynamics of the sintered 
rare-earth ferromagnet Nd 2 Fei 4 B in a transverse field, a room-temperature realization of 
the Random Field Ising Model. The power-law exponents and scaling forms of Barkhausen 
noise events as a function of temperature and transverse field demonstrate that the system 
can be placed into two disjoint regimes depending on tunable disorder: a strong-randomness 
regime in which the random-field-induced pinning dominates and a thermal regime in which 
fluctuations dominate. In the strong randomness regime, the critical exponents are close 
to mean-field predictions for heavily disordered systems. Moreover, the deep free-energy 
minima allow for domain-wall oscillations for extended periods of time. Having demon¬ 
strated the feasibility of using Barkhausen techniques to access the domain dynamics in 
this bulk random-field magnet, it should in principle be possible to extend this technique 
to single-spin, dipole-coupled, random-field magnets like LiHo 3 ;Yi_j.F 4 and Mni 2 -ac, where 
the random-field dynamics are enriched by the addition of quantum tunneling terms to the 
Hamiltonian. 
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